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Abstract 

We derive a two-band finite-dimensional model for description of the condensate tunneling in an 
accelerating optical lattice, taking into account the fine Bloch band structure. The model reveals 
a very strong dependence of the final band populations on the initial populations and phases. 
Most importantly, additionally to the known asymmetric dependence on the nonlinearity, there is 
also a notable asymmetry in the sensitivity of the tunneling probability to the nonliearity-induced 
initial population of the Bloch band to which the tunneling takes place. This fact can explain 
the experimentally observed unexpected independence of the upper-to-lower tunneling probablity 
on the nonlinearity. Finally, we compare the predictions of the two-band model with that of the 
well-known nonlinear Landau-Zener model and find disagreement when the two bands are initially 
populated. The disagreement can be qualitative and reveals itself even for a negligible nonlinearity. 
However, the two models agree remarkably well if just one band is populated initially. 

PACS numbers: 03.75.Lm, 03.65.Xp, 03.75Kk 
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I. INTRODUCTION 



Since the successful realization of Bose-Einstein condensation (BEC) in an optical lattice 

n q 

|1|, |2|] the experimental efforts are directed to explore and understand the behavior of BEC 
in a periodic potential. Many of the BEC phenomena in an optical lattice have an analog in 
the solid state physics. For instance, Bloch oscillations and Landau-Zener tunneling of BEC 
were experimentally observed SQ^ in an accelerating lattice. Tlre atonrrc interactions 
in BEC, however, give rise to the essentially nonlinear effects such as the asymmetry of the 
tunneling probability between the two adjacent bands 5,0] and the Landau and dynamical 
instabilities . Accordingly, two essentially different regimes of BEC dynamics 

in an optical lattice can be identified p, the "instability regime", which realizes for 
a small to intermediate values of the lattice acceleration and the "Landau-Zener regime", 
which realizes for larger accelerations. When the acceleration is small, the condensate slowly 



crosses the ec, 
parameter 
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Hf the Brillouin zone allowing the unstable modes to dominate the order 
. For instance, in one of the bands the modulational instability will 
result in the bright soliton formation even in a repulsive BEC jl^j ]. 

For the "Landau-Zener regime" a reduced two-level model (the nonlinear Landau-Zener 
model) was proposed in Refs. 15, ^| to account for the effect of the atomic interactions 
on the tunneling probability of BEC (see also Refs. jl^l. Il7| for further development). This 
model is derived for a weak optical lattice by an analog of the method of free electrons in 
;he solid state theory (see for instance, Ref. |lq]) and is essentially the Landau-Zener model 
191 1 modified by the nonlinear terms due to the atomic interactions in BEC. 

Tne predictions of tire nonlinear Landau-Zener mode, were experimentally tested fl ^. 
In accordance with the theoretical prediction the atomic interactions in BEC make the 
tunneling probability asymmetric, but still a significant disagreement was discovered. Hence, 
a more accurate model is in order for description of the condensate tunneling. 

To account for the modulational instability present in one of the Bloch bands a spatio- 
temporal model was proposed recently [ll|. Essentially, it consists of two linearly coupled 
nonlinear Schrodinger equations for the slowly varying amplitudes of the two Bloch waves 
bordering the Brillouin zone edge (i.e. at k = £;#). The key role of the modulational 
instability of Bloch wave amplitudes in the transition from the reversible tunneling to the 
asymmetry in the band populations is revealed. Therefore, this model is suitable for descrip- 
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tion of the "instability regime", when the Brillouin zone edge crossing time is comparable 
to the characteristic time of development of the modulational instability. However, its ap- 
plicability is restricted to the case of a very small lattice acceleration and thus it cannot be 
applied to the "Landau-Zener regime" (as defined above). 

The aim of this paper is to propose an improved finite-dimensional system (the two-band 
model) for BEC tunneling in an accelerating optical lattice. We take into account the fine 
structure of the Bloch bands, which is of crucial importance if the initial conditions are 
arbitrary. Only the "Landau-Zener regime" is considered. Hence, we can safely neglect the 
spatial inhomogeneity of the Bloch wave amplitudes. We expand the BEC order parameter 
over the Bloch waves with a time-dependent wave index. Thus our method is similar to 

ml 



Houston's approach [2Cj to the accelerating electrons in the solid state theory. The resulting 
finite-dimensional model can also have application in the study of ultracold atoms in a 
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periodic potential submitted to a constant external force 

There are two main results. First, the two-band model reveals a dramatic dependence of 
the final Bloch band populations on the initial populations and the phase difference. Most 
importantly, additionally to the asymmetry in the dependence of the tunneling probability 
on the nonlinearity strength, we reveal also the asymmetry in the sensitivity of the tunneling 
probability to the population of the adjacent Bloch band (i.e. to which the tunneling takes 
place). This fact, combined with the nonlinearity-induced population of the adjacent band, 
can explain the independence of the upper-to-lower tunneling probability on the nonlinearity 
suggested by the experimental results of Ref . P . Second, there is a significant disagreement 
between the predictions of the two-band model and the Landau-Zener model if the two 
bands are initially populated. The disagreement can be qualitative and reveals itself even 
for a negligible nonlinearity. The two models agree in the predictions if only one Bloch band 
is populated initially. 

The paper is organized as follows. In the next section a derivation of the two-band 
model is presented. Some details of the derivation are relegated to appendix |XJ Then, in 
section ITTT] the case of a weak lattice is considered. The predictions of the two-band model 
are compared to those of the nonlinear Landau-Zener model in section IIVI The two-band 
model in the diabatic basis is given in appendix [Bj Section IS contains discussion of the 
asymmetric BEC tunneling. Finally, a general discussion and perspectives are contained in 
the concluding section fVTl 
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II. DERIVATION OF THE TWO-BAND MODEL 



We consider a quasi one- dimensional BEC, i.e. the transverse width is much smaller than 
the length of the condensate. For not too large number of atoms the transverse dynamics of 



BEC is that of a quantum particle confined to the lower energy level of the parabolic trap 24] . 



In terms of the s-wave scattering length a s , the transverse oscillator length £± = a/ h/muj± 
and the linear condensate density nm — N^/d, where is the number of atoms per 
lattice site, the condensate can be considered as quasi one-dimensional under the condition 

n 

n\D a s C 1 25]. In the experiments the usual values are: d ~ l^um and a s ~ lnm. Hence, 
to satisfy the quasi one- dimensionality condition one has to have iV d <C 1000. The BEC 
order parameter $(fj_, x, t) then can be approximated as a product of the wave-function for 
the ground state of the transverse parabolic trap and a function describing the longitudinal 
motion of the condensate: 

Hr ± , X ,t) = ^e W ! [ -±^( X ,t). (1) 

For the reduced order parameter ip(x,t) the Gross-Pitaevskii equation describing the con- 
densate in an accelerated optical lattice reads 



W = + cos(2A; L [x - x (t)]) + ^x 2 \ $ + (2) 

where g = 4nh 2 a s /m, k^ is the optical lattice index and Xo(t) describes the acceleration 
of the lattice. Here the order parameter is normalized as follows d~ x dx \ip\ 2 = 1 with 
d = n/ki being the lattice spacing. The following dimensionless variables will be used below: 

t=^-t, x = 2k L [x-x (t)}, v = -^, c = , A 2 = (2£||A; L )~ 2 , (3) 

where Er = h 2 k'j j /(2m) is the recoil energy and t\\ = ^Jh/muj\\. We have 

1 ( , dig \ ~ ( /~\ , ^ r~ i ~ 12^ 7 I 7|2 



idpp = ~ [-id* - j ip+ Uqcos(x) + y[x + x ] ) ^ + c\ifj\ 2 ip } (4) 



where x = 2kix Q and 



t / ., 
i f ( dx§ 



*P = exp{ -- J drl ] 



o 

r2vr , _ , 7 i 



The normalization now reads (2ti) 1 J Q * dx \ip\ 2 = 1. As we will use only the dimensionless 
variables of equation @, we drop all the tildes from x, t, and ip, for simplicity. 
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We assume that the quasi-momentum spread of the condensate in the optical lattice is 
negligible compared to the width of the Brillouin zone. One necessary condition is that the 
optical lattice be long enough, i.e. the parabolic trap contains many lattice periods: Ac 1. 
Then, the effect of the parabolic trap is negligible. We set A = 0. 

An additional condition is that the time of the Brillouin zone edge crossing is much less 
that the characteristic time of the modulational instability development. This sets a lower 
limit on the acceleration a. The latter time can be estimated as follows. In the slowly varying 
envelope approximation 11. Il^ the spatial Bloch amplitude will be governed by a nonlinear 



Schrodinger equation with an effective interaction coefficient cxnn- The characteristic time 
of the instability development is the inverse of the maximal growth rate tui ~ l c l _1 - 

Formation of gap solitons is also possible for a nonlinearity of the order of the gap between 
the Bloch bands [26|. In the following, however, we consider only the weak nonlinearity limit 
c<!) , used in the experiment [a]. 

Under the conditions that the condensate has a narrow initial momentum spread and 
the time of the Brillouin zone edge crossing is much less that the characteristic time of the 
modulational instability development the solution to equation (J3J) can be expanded in the 
basis of the Bloch waves with a definite (time-dependent) index q(t) and an overall linear 
phase due to the acceleration: 



n=l 

where v(t) = xo and ^ n>q (x) = e %qx u n) q{x) satisfies 

|-^+VoCOs(ar)J*„ )ff (a7) = E n (g)* n ,,(ar). (6) 
Due to the normalization and orthogonality of the Bloch functions, 

7T OO 

(27r) _1 f dx^ k (x)^ m} k(x) = S nm , the amplitudes satisfy the condition \A n \ 2 = 1. 

— 7r n=l 

It is shown below that, in a vicinity of the Brillouin zone edge q = ks, for a small 
acceleration and under the condition E 3 (kB) — E 2 (kB) 3> E 2 {kB) — Ei(ks) one can neglect 
the effect of the upper bands on the tunneling between the first two Bloch bands. The 
nonlinear term with c <C v o can be taken into account by using a perturbation theory. 

The expansion (0) is in the spirit of Houston's approach in the theory of electrons accel- 
erating in a crystal lattice j^. It takes into account, for instance, that in the experiments 



the condensate was in a Bloch state before the acceleration was turned on. 



A. The linear limit c = 



For the moment let us drop the nonlinear term and consider the linear limit of equation 
(@J). The contribution from the nonlinear term is accounted later on. Substituting the 
expression (J3J) into equation (JU) and dropping the overall imaginary exponent gives 



dq dv 
dt + dt 



X ^ A n U n:q (x) + ^ i ~^ !Lu n,q(x) + Z^j ^ A n d q U n , q (x) 



= y^ j E n (q)A n u ri!q (x). (7) 

Setting g(t) = q — v(t) to cancel the non-periodic in x term and projecting the resulting 
equation onto the basis function u n>q (x) we get an equation for the amplitude 

E n (q)A n + ia(t)^T A m (u n , q \d q u m , q ) , (8) 



dt 



where a = x . Here the inner product over the lattice period is used: (/1I/2) = 
(27r) _1 f* n dx fZ(x)f 2 (x) with {u n , q \u m>q ) = S n>m . 

In appendix 1X1 it is shown that for any periodic potential v = v(x), with v(x) being an 
even function, u*(x) = u n ^ q (—x) and, hence, the diagonal inner product (u ntq \d q u niq ) = 0, 
while the non-diagonal is given as 

K (0) = (U \3U ) = -j ( Un >l\ dxUm ^ = -j ( U n,q\-jr\Um,q) , (Q) 

KnmW - {U nM U m , q) I ^ _ ^ , ^ _ ^ )]2 , U ? m. W 

The coupling coefficient K nm is real and satisfies n nm (<l) = —^mn( ( l)- 

Note that the coupling coefficient between the Bloch bands is inversely proportional to 
the band level spacing E m {k) — E n (k), since n nm (k) = O (vo/[E m (k) — E n {k)] 2 ). Thus, for 
a small acceleration and an optical lattice satisfying E 3 (ks) — E 2 (ks) E 2 (kB) — -&l(/cb), 
the influence of the upper bands on BEC tunneling between the first two Bloch bands at 
the Brillouin zone edge is negligible. For instance, the above condition is satisfied by a weak 
optical lattice (see section ITTT|) . 

Dropping the terms corresponding to the upper bands from equation (JSJ) and making the 
transformation A\ = exp j— i dr E(q(r))\ a\ and A 2 = exp j— i f* dr E(q(r))\ (— 1)0,2, 
with E = (Ei + E 2 )/2, we obtain the two-band model for the linear case: 

i ~dT = ~ £ ^ ai + aKl2 (^) a2 ' ( 10 ) 
6 



i^j- = e(q)a 2 + aKi 2 (q)ai. (11) 

Here we have used that K2i(<?) = — «i2(?) and set e(q) = (E 2 (q) — E\(q))/2. Recall that 
a(t) = Xo is the lattice acceleration and q(t) = qo — xq is a time-dependent parameter. 
The system of equations (|lUj) - (jllj) is consistent with the two-band approximation, since it 
conserves the total number of atoms in the two bands: |ai(i) | 2 + \a,2(t) \ 2 = 1. 

The linear two-band model ([10 )1 - ([11 )1 is important in its own right, since it can have direct 
application in the study of ultracold atoms in a periodic potential submitted to a constant 



external force 
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B. Weak nonlinearity: c <C v q 



In the experiments on BEC tunneling in an optical lattice ja, |y] the nonlinearity satisfies 
the condition c <C Vq. In this case, the nonlinear term can be treated as a perturbation of 
the optical lattice potential. Here it is pertinent to mention that the nonlinear term was 
also treated by introduction of an effective optical lattice potential 0,0]. This treatment 
is of limited validity. 

Under the condition c <C v the nonlinear term preserves the qualitative Bloch band struc- 
ture, but modifies the Bloch functions and the band energies. The nonlinear term serves as 
an additional lattice potential (in general, time- dependent), the same for all eigenfunctions. 
In fact, if we define the nonlinear modes u n) k(x,t) as follows 



■(—id x — k) 2 + vq cos(x) + c 



^ ] ^m(t)u mi k(x, t) 



m=l 



u n ,k{x, t) = E n (k, t)u n ,k(x, t), (12) 



then the effect of the nonlinear term can be treated by using the usual perturbation theory 
for the eigenvalues E n (k,t) of a perturbed ((A;, independent) linear operator. Notice that 
the nonlinear modes u n>q are orthogonal, (u mt k\u n ,k) = 5n,m, and for c = coincide with 
the respective Bloch functions. Thus, for small c, they constitute a complete orthogonal 



basis (in this respect we generalize the approach of Refs. 
perturbation theory we get 
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). In the first order of the 



E n (k,t) = E n (k) + c(u* (x) 



A ™(fy 



m=l 



u n>g (x)) + 0(c 2 ), 
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with (...) = (27r) _1 dx(...), while the eigenfunctions u njq (x) = u n>g (x) + 0(c) which is 
a sufficient approximation for the following. In the two-band approximation we get the 
following expressions for the energies of the first two Bloch bands (the variable t is omitted 
for simplicity): 

Ei(k) = E 1 (k) + c { X ii(k)\M 2 + Xn(k)\M 2 + ovi{k)A x J% + a{ 2 (k)A\A 2 } , 



E 2 (k) = E 2 (k) + c { X 22(k)\M 2 + Xi2(k)\M 2 + <T2i{k)A 2 A\ + a^k)^} . (13) 

Here 

Xnn = (|Mn,fe| 4 ), Xl2 = ( |«l,fc | 2 \u 2 ,k \ 2 ) , 0"l2 = (|«l,fc| Vl,fc«2,fc) , °"21 = ( | ^2,fc 1 2 «2,fc^ ife ) . 

The terms with er nm can be neglected, since at the Brillouin zone edge these coefficients are 
exactly zero due the properties of the Bloch waves. Numerics shows that they grow linearly 
in k — ks- 

The nonlinear generalization of the system ()10 j) -()ll j) is derived by using the nonlinear 
Bloch waves ^ n ,k = G lkx u n ^ = *& n ,k + C( c ) i n the expansion (j5J). For a small acceleration 
the term proportional to ac can be neglected, thus the inter-band coupling coefficient K12 
in the nonlinear system is given by the same formula as in the linear case. Setting now 
A\ = exp |— i J ' dr F(g(r))| a\ and A 2 = exp |— i J * dr F(g(r))| (— ia 2 ), with F = {E x + 
E 2 )/2 + CX12, we obtain the nonlinear two-band model: 

= {~ £(yq " > + C7ll ^)l ai ' 2 } fll + aK ^(q)a 2 , (14) 

= {s{q) + C7 22 (g)|a 2 | 2 } a 2 + otK 12 {q)ax. (15) 
Here we have used the modified Bloch energies from equation (|13|) and defined 

Inn = Xnn ~ Xl2 = (\Un,k\ 4 ~ \ U l,k\ 2 \ U 2,k\ 2 ) 1 U = 1,2. (16) 

The system fll4|) - (fl3j) depends on the fine /c-structure of the Bloch bands, since q(t) scans 
through the edge of the Brillouin zone. Hence, the expressions for the Bloch waves in a 
vicinity of the Brillouin zone edge have to be available. For a general optical lattice the 
coefficients e(q), ni 2 (q), and 7„„(g) can be computed numerically. However, such expres- 
sions are available in the analytical form (see below) for a weak optical lattice used in the 
experiments 0, 6|. 



III. WEAK LATTICE APPROXIMATION 

Similar to the free electrons approach in the solid state theory ^| , in the vicinity of the 
Brillouin zone edge the Bloch waves from the two lowest bands of a weak optical lattice can 
be sought in the form 

$ = ae iqx + be i(q - 2kB)x . (17) 
For the cosine lattice v(x) = f cos(x) we get 

ExM = q l + k B (k B -q) T e(q), e(q) = (k 2 B {k B - qf + |^ . (18) 

The corresponding Bloch waves read 

m) e e e + i ^ qje 

where v{q) = 2[k B (k B — q)+e(q)]/v . Then the remaining two coefficients of the two-band 
model are 

. 1 v(q) Av 2 (q) 

KlM = Mv)TT^Wy 7(g) " 7 - (9)= [i + , 2 (g)p - (20) 




0.25 0.5 0.75 1 0.25 0.5 0.75 1 



q q 

FIG. 1: The Bloch energy of the first two bands, panel (a), and the inter-band coupling coefficient 
K12, panel (b), vs. the band index q. Numerically found values are given by the solid lines and the 
analytical weak-lattice approximations are given by the dashed lines. Here vq = 0.1375. 

For the linear case, the analytical formulae in equations (|18|) and (|20|) are in good agree- 
ment with the exact numerical results for the experimentally used values j^lly (we use 
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vo = 0.1375). In fig. [I] the analytical energy levels Ei^{q) (left panel) and the inter-band 
coupling coefficient Ku{q) (right panel) vs. their numerical values are given. The coupling 
coefficients to the higher Bloch bands are small as compared to k 12 , in accordance with 
the theory of section 111 Al For instance, we have for the maximum value: n^ks) < 0.2, 
n = 1,2, for v = 0.1375. 

1.2 
1 

0.8 
jj 0.6 
0.4 
0.2 

0.2 0.4 0.6 0.8 1 

q 

FIG. 2: The nonlinear coefficients 711 (?) and 722(9)- The numerically computed coefficients are 
given by the solid lines and the analytical weak-lattice approximation is given by the dotted line. 
The dashed lines are the scaled analytical coefficients as indicated in the text. Here vq = 0.1375. 




The weak lattice approximation cannot, however, capture the weak asymmetry in the 
two nonlinear coefficients 711 (?) and 722(5), since their analytical values always coincide. In 
fact, in the vicinity of the Brillouin zone edge, the analytical approximation (dotted line in 
fig. |2J) is close to the algebraic average of the exact numerical coefficients. To have a better 
approximation the analytical coefficients must be scaled by the maximum value (at q = ks)- 
For instance, 711 (?) = 1.1327(?) and 722(9) = 0.8427(?) for t> = 0.1375 (the dashed lines in 

fig. ED. 

IV. COMPARISON WITH THE LANDAU-ZENER MODEL 

To make a quantitative comparison with the nonlinear Landau-Zener model 

.da [at c n 2 |2\1 , v o-, 
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.d& f at i c /l; |2 I i2\l 7 , 



nrf 

introduced in Refs. [la, I16j . one has to relate the amplitudes a and 6 of the incident and 
Bragg-scattered waves with the amplitudes of the Bloch waves from the first two bands. 
There are two approaches. The first one is based on the fact that the adiabatic energy levels 
of the quantum-mechanical Hamiltonian 

(at vq \ 
2 2 ) 

associated with the system (j21j) -(|22 j) . mimic the structure of the first two Bloch bands. 
Therefore, the Bloch amplitudes are obtained by projecting the vector (a, b) onto the in- 
stantaneous adiabatic eigenvectors of H jsjj. We obtain the following relation 

xa — b a + xb 

ai = /i~t — 2 ' ° 2 = per — 2 ' (24) 

V 1 + M V 1 + X 1 

where k — {E — at)/v with E = a/ a 2 t 2 + Vq. The amplitudes a\ and 02 in equation 
correspond to the instantaneous adiabatic energy levels E\p = ^fE/2, thus they can be 
taken as the occupation amplitudes of the first two Bloch bands according to the nonlinear 
Landau-Zener model (however, one can get only the absolute values of the amplitudes 0^2, 
their phases are left unknown due to arbitrary phases of the eigenvectors). 

The second approach is to use the analytical results for a weak lattice and equate directly 
the order parameters corresponding to the two models. The two-band model p5j1 - (jlfijl 
corresponds to the order parameter 

V> = e^+^W* [ ai (t)* liq(t) (x) - ia 2 (t)* 2 , q(t) (x)] , (25) 

with 8(t) being an unimportant phase and q(t) = qo — v(t). The nonlinear Landau-Zener 
model ()21jl - ()22jl corresponds to the following order parameter 

= e~ 3ict/2 [a(t)e ikx + b(t)e i{k - 2kB)x ] . (26) 

It is straightforward to relate the band amplitudes using the expression for the Bloch waves 
([12)) and associating k = q . The common x-independent phase can be discarded. We have 







: 





— 1 —iv 
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)l 









(27) 



The transformation given by equation ()27|) is linear and unitary. Asymptotically it trivializes. 
Indeed, v — ► oo as t — ► — oo, while v — > as t — ► oo (here a > 0). We get 

a(— oo) = ai(— oo), 6(— oo) = — ia 2 (— oo), a(oo) = — za 2 (oo), 6(00) = — ai(oo). (28) 

1 

0.75 
CL 0.5 
0.25 


0.1 0.2 0.3 0.4 0.5 

a 

FIG. 3: The tunneling probability vs. the lattice acceleration according to the Landau-Zener 
formula (solid line) and in the two-band system (points) in the case of negligible nonlinearity and 
trivial initial condition (i.e. just one band being populated initially). Here vq = 0.134 and c = 0. 

The overall conclusion to be argued below is that, whatever the way one relates the 
two models, the predictions of the nonlinear Landau-Zener model (|2i p -([22| ) and the two- 
band model (|15 p -()16| ) disagree significantly in general, if the initial condition is non-trivial, 
i.e. both amplitudes (a and b and, respectively, a\ and a<i) are non-zero initially. The 
disagreement is apparent even for a negligible nonlinearity, which case is treated below in 
detail. A remarkable agreement between the two models if just one band is populated 
initially (i.e. either a\ or a 2 is zero initially) must also be stressed. 

When initially the atoms populated just one band the two models give the same results. 
For instance, for c = the probability of tunneling, defined as P = |a n (oo)| 2 for a n (— 00) = 0, 
is given by the Landau-Zener formula P = e~ nv o^ 2a " > |l9( which is illustrated in fig. |H1 

When both bands are populated initially, the two models do not agree in the prediction 
of the final populations. This is illustrated in fig. EJ where we plot the time-dependence 
of the the second band population for different initial conditions but with the majority of 
atoms populating initially the first band. Note that in panel (a) the two models give the 
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FIG. 4: The population of the second Bloch band vs. time in the two-band model (solid lines) 
and in the nonlinear Landau-Zener model (dashed lines). We have used vq = 0.1375, a = 0.036, 
and c = and the following initial conditions: (a) (01,02) = (1>0), (b) (01,02) = (V0.99, VO.01), 
(c) (oi,o 2 ) = (VbT95,\/0T05), and (d) (oi,o 2 ) = (V095, e~ in / 5 V0W>)- The initial amplitudes a 
and b for the Landau-Zener model and the corresponding functions ai(t) and 02 (t) were computed 
using equation (|27|) . The initial time was to = —200. 



same results if related by equation (j2Zj). However, the same equation leads to significant 
disagreement in panels (b) and (c). Similar conclusion can be drawn using the first relation 
scheme for the two models, i.e. equation plj) . 

Two more conclusions are apparent from fig. 0] When the first band is initially domi- 
nantly populated, a weak population of the second band makes a strong effect on its final 
population, (compare panels (a) and (c)). Moreover, at the fixed initial populations, the 
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initial phase difference can have a significant effect on the final populations (compare panels 
(c) and (d)). The same conclusions are true for the first band. 




-50 -25 25 50 
t 

FIG. 5: The populations of the Bloch bands vs. time in the two-band model (thick lines) and in 
the Landau-Zener model (thin lines). The solid lines give |ai| 2 . The initial time was to = —200. 
Here vq = 0.1375, a = 0.036, and c = 0. We have used the following initial condition (01,02) = 
(i\/0.15, V0.85). The initial amplitudes a and b for the Landau-Zener model and the corresponding 
functions a\(t) and 02 (i) were computed using equation (|2*7|) . 

The significant disagreement between the two models can be manifested by using a larger 
initial population of the less populated Bloch band. In this case the disagreement is quali- 
tative, see fig. El 

The transformation ^Tf\ can be used to derive the two-band model in the diabatic basis, 
i.e. in terms of the amplitudes of the incident (a) and Bragg scattered (b) waves. The 
resulting system is given in appendix [Bl The main feature of the two-band model in the 
diabatic basis is that both the sweep function and the coupling coefficient are complicated 
functions of the band index q with the coupling coefficient being complex valued. This 
latter fact can explain the disagreement in the predictions of the two-band model and the 
nonlinear Landau-Zener model when the two bands are initially populated. 
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V. ASYMMETRIC TUNNELING FOR c / 



The tunneling probability of BEC at the edge of the Brillouin zone in an accelerated 
optical lattice was found to have an asymmetric dependence on the nonlinearity coefficient 
c l5j. The asymmetry in the tunneling probability was first theoretically predicted in Ref. 
161 ] within the nonlinear Landau-Zener model. The two-band model ()15J) - (fT()|) confirms the 
existence of the asymmetric tunneling, see fig. |3 



0.55 
0.5 
CL0.45 
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FIG. 6: The tunneling probabilities in the two-band model (solid lines) when just one band is 
initially populated. The upper line gives the lower-to-upper tunneling probability and the lower 
one gives the upper-to-lower tunneling probability. For comparison, the tunneling probabilities in 
the nonlinear Landau-Zener model are also given (dashed lines). Here vq = 0.1375 and a = 0.036. 

r-i 

The experimental results (see fig. 2 in Ref. p) also suggest that the upper-to-lower tun- 
neling probability, in contrast to the lower-to-upper one, is unaffected by the nonlinearity 
(there is, however, a large error due to imperfections in the experiment). Below we argue 
that, indeed, there is an additional asymmetry in the sensitivity of the tunneling probability 
to the initial nonlinearity- induced population of the adjacent Bloch band, which, together 
with the asymmetric dependence on the nonlinearity coefficient, can give a complete expla- 
nation of the experimental results. 

In section HVl (see fig. we have found that the final populations of the Bloch bands are 
affected by the initial band populations and phases. Due to imperfections in the experimental 
preparation of the initial state, there is always a small population of BEC atoms in the 
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adjacent Bloch band. The condensate is initially prepared at the center of the Brillouin 
zone (k = 0) with a large fraction of atoms populating one of the two Bloch bands. Besides 
the experimental imperfections, nonlinearity (i.e. the atomic interactions in BEC) also 
contributes to a small number of BEC atoms in the adjacent band, since the nonlinearity 
couples the Bloch bands in the higher orders of the perturbation theory discarded in section 

urn 

To find out the effect of the nonlinearity- induced population of the Bloch band on the 
tunneling probability to this band, we choose to use the following initial condition aj = y/rc, 
where j = 2 for the lower-to-upper tunneling and j = 1 for the upper-to-lower tunneling 
(we do not know the experimental initial share |aj| 2 of the number of atoms in the adjacent 
Bloch band, however, this number grows with the nonlinearity coefficient c [31]). Since the 
Bloch waves are real for k = the initial phase difference between the amplitudes a\ and 
ci2 of the two-band model is close to zero and does not affect the general picture. For the 
coefficient r we have used the following set of values {0.01; 0.05; 0.075; 0.1; 0.15; 0.2}. The 
results of the numerical simulations of the two-band model are presented in fig. The 
probability of tunneling (defined here as the final population of the Bloch band to which the 
tunneling takes place, neglecting a very small initial population of this band) is given vs. 
the nonlinearity coefficient c. The asymmetry in the spread of the numerical curves with 
different r in the two cases of tunneling is apparent. Here we should mention that fig. [7|can 
provide only the qualitative account of the asymmetry in the sensitivity of tunneling to the 
initial Bloch band populations, since the weak lattice approximation developed in section 
IIHI applies only in the vicinity of the Brillouin zone edge. 

It is easy to explain the tunneling probability change if the Bloch band to which the 
tunneling takes place is populated with a small number of atoms. Indeed, the nonlinear 
terms in the two-band system (|14 p -()15| ) effectively raise the band energy levels. Therefore, 
in the case of the lower-to-upper tunneling a small population of the adjacent band increases 
the energy gap, while in the case of the upper-to-lower tunneling it decreases the energy gap, 
thus the tunneling probability is decreased in the first case and increased in the second one. 
However, this argument cannot explain the fact that the strength of the sensitivity to the 
population is quite different in the two cases. 

In conclusion, a very small initial population of the Bloch band to which the tunnel- 
ing takes place (in fact for r = 0.2 we get \aj\ 2 = 0.12) has a dramatic effect on its final 
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population in the case of the upper-to-lower tunneling as compared to that in the case of 
the lower-to-upper tunneling (if one consider the relative change in the tunneling proba- 
bility when increasing r). This fact, together with the known asymmetry of the tunneling 
orobability, provides a complete explanation of the experimental results presented in Ref. 
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FIG. 7: The lower-to- upper (a) and upper-to-lower tunneling (b) in the two-band model with a 
weak nonlinearity- induced initial population of the adjacent Bloch band (the upper band in case of 
(a) and the lower one in case of (b)). The probability of tunneling (the final population of the Bloch 
band to which the tunneling takes place) is given vs. the nonlinearity coefficient c. The initial 
conditions are (a) [01,02] = [VI — rc, y/rc\ and (b) [01,02] = [y/rc, VI — fc\. The dashed lines give 
the tunneling probablity for r = 0. The solid lines have r = {0.01; 0.05; 0.075; 0.1; 0.15; 0.2} (the 
solid lines deviate from the dashed one as r increases). Here vq = 0.134 and a = 0.036. The time 
interval is [—60,60]. 



VI. CONCLUSION 

In the present work a finite-dimensional model has been derived for the condensate tun- 
neling at the edge of the Brillouin zone in an accelerating optical lattice. The fine Bloch 
band structure was taken into account. A special attention was paid to the case of a weak 
lattice, which could have seemed a redundant task due to the existence of the nonlinear 
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15. 161. However, it turns out that 



Landau-Zener model proposed for this case previously 
the two-band model derived in this paper and the Landau-Zener model disagree significantly 
even in the case of a negligible nonlinearity if both adjacent bands are populated initially. 
Only in the special case of just one Bloch band being initially populated the two models 
agree in the prediction of the final populations. In particular, the Landau-Zener result 3| 
is reproduced by our model. 

The two-band model confirms the asymmetry in the tunneling probability dependence on 
the nonlinearity, predicted theoretically in Ref. |16] and discovered experimentally in Ref. 
0|. However, besides the qualitative agreement between the theory and the experiment, a 
quantitative disagreement was also revealed. For instance, the experimental upper-to-lower 
tunneling seems to be unaffected by the nonlinearity |5|]. We have proposed an explanation 
for this fact based on the discovered strong dependence of the final Bloch band populations 
on the initial populations. Indeed, additionally to the asymmetry in the dependence of 
the tunneling probability on the nonlinearity strength, we find also the asymmetry in the 
sensitivity of the tunneling probability to the initial populations: the upper-to-lower tun- 
neling probability is stronger affected by a variation of the initial Bloch band populations 
than the lower-to-upper tunneling probability. For instance, about one per cent of BEC 
atoms populating the lower Bloch band due to nonlinearity-induced Bloch band coupling 
significantly modifies the upper-to-lower tunneling probability, whereas the lower-to-upper 
tunneling probability is only weakly affected by the same initial population of the upper 
Bloch band. 

Though it is in the dependence of the final Bloch band populations on the initial ones 
where our model disagrees with the nonlinear Landau-Zener model, nevertheless, both mod- 
els show strong dependence of the final Bloch band populations on the initial populations. 
Thus an experimental study and/or numerical simulations of the full GPE equation with 
various initial Bloch band populations is in order to confirm the strong dependence of the 
final populations on the initial conditions and compare with the two models. 

So far we have discussed the so-called "Landau-Zener regime" of BEC tunneling neglecting 
the spatial modulation of the Bloch wave amplitudes. Recently the modulational instability 
of a Bloch wave amplitude was found to result in the irreversibility of tunneling and the 



asymmetry in the final band populations 



1 If ] . A unified approach, capable to cover the 



two identified regimes, the "Landau-Zener regime" and the "instability regime" can also be 
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developed. The solution of this problem is left for a future publication. 
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APPENDIX A: THE COUPLING COEFFICIENT BETWEEN THE BLOCH 
BANDS 

The inner product (u n<q \d q u mtq ) can be easily evaluated by using the properties of the 
Bloch functions u n ,q(x). For generality, consider an arbitrary periodic potential v = v(x). 
Define the Hermitian operator C q = \{—id x + q) 2 + v(x). We have then £ q u n , q = E n (q)u njq 
and, hence 

CqdqU ntq = E n (q)d q U n) q + ~^~ U n,q + ^xU n ,q ~ qU n ,q- (Al) 

Using this result and orthogonality of the functions u n>q we get 

E n {q){u nA \d q U mA ) = (u n: q\Cq\dqU mt q) = E m {q){u n>q \dqU m>q ) + i(«n,g \d x U m , q ) . 

Therefore, for m ^ n we obtain 

(u n ,q\d qUm ,q) = _<^^-2^S_. ( A 2) 

Now, let us evaluate the product on the r.h.s. of equation (jA2|) . We have 

r a - w ( \a - dv ^ 

L, q O x U mq ^- , m\q )O x Um,q ^m,q- 



Hence, 



■En\,Q) (u n q \d x Um,q) (^n,q\C'q\d x Urn,q) E m (q^J (u n,q\d x U m q ) (Un,q\ |^m,<j) 



with the result 



K nm - {U nt q\d q U mt q) - . _ - _ y , U * m. { A6) 
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What is left is to consider the diagonal inner product (u n!q \d q u njq ) . Here we use that the 
potential is even function: v(—x) = v(x). Then the functions u n , q (x) can be selected to 
satisfy u^ q (x) = u ntq (—x). We get 

(u ntq \d q u n:q ) = ^- J dx u n , q (-x)d q u n , q (x) = J dx (-d q u n>q (-x)) u n>q (x) 

— TV —7T 
7T 

= ^ J dx (-d q u n>q (x)) u nA (-x) = -{u nA \d q u n , q ), 

— IT 

where to replace the g-derivative we have used the independence of the inner product on q: 
(u n ,q\un,q) = 1- Hence (u n ^ q \d q u n>q ) = 0. A similar reasoning for (u n , q \(— i)d x u m>q ) can be 
used also to establish that n nm (q) is real and that K nm (q) = —K mn (q). 

APPENDIX B: THE TWO-BAND MODEL IN THE DIABATIC BASIS 

Unitary transformation (J27)) for the two-band model can be interpreted as a relation 
between the adiabatic basis (01,02), where the amplitudes are the band populations, and 
the diabatic basis (a, b), where the variables have the physical meaning of the amplitudes of 
the incident and Bragg scattered waves. For simplicity, we consider below only the linear 
case, c = (the nonlinearity does not pose any problems, but complicates the presentation). 
After simple calculations one arrives at the following system in the diabatic basis: 

i^ = E(g)o + K(g)6, (Bl) 

i^ = -E{q)b + lC(q)a, (B2) 
where the coefficients are as follows 

l + v 2 (q) l + v 2 (q) e(q) 1 + v 1 {q) 

Here the functions e(q) and v{q) are taken from section ITTT1 (The system (|B1[) - ([B3[) was also 
verified numerically against the two-band model.) 

It is easy to see that E(g) is an odd function of kg—q, while K(g) is an even function of this 
variable. The following asymptotic properties of the coefficients can be easily established. 
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Recalling that q = — a we get as t — > ±00: 

E(g)--y, K(g)-|. (B4) 

At the Brillouin zone edge we get E(/cb) = and K(&b) = Vq/2 + ia/vo- 

However, an important fact is observed about the system (jBl|) - (jB2|) : the coupling coef- 
ficient K(g) is a complex valued function. Therefore, to obtain an equivalent Landau-Zener 
model (with a real coupling coefficient), one has to perform the following time- dependent 
phase transformation 



a = a 



(LZ) e »e(,) j b = & (LZ) e - i e( 9 ) ) e(g) = arg{K(g)}/2. (B5) 



For the phase-transformed amplitudes we obtain: 

da( LZ ) 



dt 

d& (LZ) 



l- 



E (LZ)^ a (LZ) + | K ( 9 )| 6 (LZ) j (gg) 



. E (LZ) (g)6 (LZ) + |K(g)|a (LZ) ) (B7) 



dt 

where E< LZ )(g) = E(q) + a&(q) (& = dQ/dq). 

It is easy to see that E^ LZ ) — > —at/2 as ±00 and E( lz )(/cb) = 0, i.e. we have all the 
necessary prerequisites for application of the arguments due to Landau and Zener to 
the system ()B6j) - (jB7)) . This explains the agreement between the Landau-Zener model and 
the two-band model when just one band is populated initially. 
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